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Abstract 

We take charged anisotropic fluid cylinder when there is no ex- 
ternal pressure acting on the fluid. This is a cylindrical version of 
the Krori and Barua's method to explore the field equations with 
anisotropic fluid. We discuss models with positive matter density and 
pressure that satisfy all the energy and stability conditions. It is found 
that charge does not vanish at the center of the cylinder. The equilib- 
rium condition as well as physical conditions are discussed. Further, 
we highlight the connection between our solutions and the charged 
strange quark stars as well as with dark matter including charged 
massive particles. The graphical analysis of the matter variables ver- 
sus charge is given which indicates a physically reasonable matter 
distribution. 
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1 Introduction 

The study of general relativistic charged compact objects is of fundamental 
importance in astrophysics. Strong magnetic fields, different kinds of phase 
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transitions and solid stellar core cause anisotropy in the fluids. However, 
charged fluids with anisotropy complicates the solution of the field equa- 
tions. Equations of state (EoS) has important consequences in such situa- 
tions. Many exact solutions have been obtained [1] by using a simple form 
of the energy-momentum tensor and assuming some symmetries. 

There have been pervious discussions of a similar nature by Evan [2J, 
Bronnikov [3] , Latelier and Tobensky [3] , and Kramer [5] . They used various 
EoS that could be written in the form p = for specific positive values 
of 7, as well as energy conservation. Some work has been done on charged 
anisotropic static matter by using spherically symmetric stars with the linear, 
nonlinear and Chaplygin gas EoS. Ivanov [6J showed that the field equations 
can be simplified by using linear EoS for a charged perfect fluid but with non- 
integrable equations. Sharma and Maharaj [7] explored the field equations 
for static spherically symmetric uncharged anisotropic fluid with combined 
linear EoS and a particular mass function. 

Charged anisotropic fluids have been discussed in General Relativity since 
the pioneering work of Bonner [8]. Ray et al. [H] investigated charged 
anisotropic spheres with Chaplygin gas EoS. Thirukkanesh and Maharaj [TU] 
generated models for charged anisotropic spherically symmetric stars by us- 
ing linear EoS as well as choosing one of the metric functions and electric field 
intensity. Horvat et al. [11] studied gravastars for charged anisotropic fluid. 
Recently, Victor et al. [12] explored solutions for the charged anisotropic 
spheres with linear or nonlinear EoS. 

Over the years, many authors have proposed various formulations to solve 
the field equations for cylindrically symmetric spacetime. Nilsson et al. [13] 
investigated cylindrically symmetric perfect fluid models. One of the authors 
(MS) [H] explored perfect fluid, static cylindrically symmetric solutions of 
the field equations by using different EoS. Sharif and Fatima [H] worked for 
the charged anisotropic cylinder but they discussed gravitational collapse. 
Som [16] explored the charged dust cylinder. However, there has been a 
little progress towards investigating charged anisotropic static cylindrically 
symmetric solutions with or without using an EoS. 

In a recent paper [T7], we have explored exact solutions of the field equa- 
tions for the charged anisotropic static cylindrically symmetric spacetime 
using Thirukkanesh and Maharaj [10J approach. Here we extend this study 
for the charged anisotropic static cylindrically symmetric spacetime by using 
Victor et al. [12] procedure. A system of differential equations for matter 
as well as electric field intensity and anisotropic pressures are solved on the 
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basis of linear and nonlinear EoS. Numerical factors depending on matching 
conditions are used with each EoS which provide relationship among charge 
distribution, pressure anisotropy and EoS. 

The outline of the paper is as follows: In the next section, we write down 
the Einstein-Maxwell field equations for the static cylindrically symmetric 
spacetime and also express this system of equations with original Krori and 
Barua's [T8] assumptions. We apply the central and boundary conditions 
on electric field intensity and radial pressure respectively to analyze these 
field equations at center and on the boundary of the cylinder. Section 3 
investigates models for linear, nonlinear and Chaplygin gas EoS by taking 
positive matter densities and pressures corresponding to the relevant EoS. 
In section 4, we match smoothly the interior and exterior metrics and bring 
adimesionality in the three models. Section 5 provides some physical features 
of these models. In particular, we discuss the stability conditions, energy 
conditions, the ven der Waals EoS [19] and the equilibrium conditions for our 
models. The last section 6 contains concluding remarks about the results. 

2 The Field Equations 

We take the static cylindrically symmetric spacetime given by [20] 

where v and \i are functions of r. The transformation dcfi = e u d6 leads the 
above equation to the following form [121 [T5] 

ds 2 = e 2 »dt 2 - e 2 »- 2v dr 2 - r 2 d6 2 - e 2 ^dz 2 . (2) 

The field equations for the charged anisotropic source are 

Ra b -\ga b R = <T^+T^), (3) 

where T^ 1 and T^™^ 1 are the energy-momentum tensors for anisotropic mat- 
ter and electromagnetic field respectively. The energy-momentum tensor for 
anisotropic fluid is 

T ir } = (P + Vr)u a U b - p t g a b + (Pt - Pr)VaVb (4) 
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satisfying u a u a = —r] a 7] a = 1, where p is the charge density, p r is the ra- 
dial pressure, p t is the tangential pressure, u a = e u 5® is the 4-velocity and 
rj a = —e^~ u 5 a is the 4-unit vector. The energy-momentum tensor for the 
electromagnetic field is 

T% m) = ^(-9 cd F ac F bd + \g ab F cd F cd ), (5) 

where F ab = A bya — A a j> is the Maxwell field tensor and A a is the 4-potential. 
The Maxwell field equations are given by 

jab 



[V=gF ab ], b = Anr^, F [abA = 0, (6) 

where J a = au a is the 4-current of the fluid element and a is the proper 
charge density. 

The field equations for the line element ([2]) become 



e 



2u-2fi 



e 2v - 2 »{v" - n") = 87tp + E 2 , (7) 
rv' 2 + rv'p' + //) = 8np r -E 2 , (8) 

e 2»-2^ u , 2 + ^ = 8npt + E 2^ (g) 

a = (reT v E)', (10) 

where prime denotes differentiation with respect to r and E = 2y/ne~ ll ^A 
stands for em part. In the system of equations (!?])- (jTOl) . there are seven 
unknowns, so we make physically reasonable choices for any two of the un- 
knowns. We take the gravitational potential e 2 ^~ 2v and the electric field 
intensity E as [TTJ 

e 2 — = 1 + (Cl ' C2)r , (11) 

1 + C\T 

2 k(3 + Cl r) 

E ~ (TT^p (12) 

where ci, c 2 and k are constants. Substituting Eqs. ffTTj) and f|T2|) in Eq.f rroj) . 
we obtain 

a^. (13) 
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This shows that there is a singularity in the charge distribution at r = 0. 
However, this choice keeps the charge distribution regular at the centre of 
the cylinder as E(r) remains finite there (1121) . 

The singularity free models for charged anisotropic static cylinder are 
constructed by taking [18J 



fi = Ar 2 , u = Br 2 + C, (14) 

where A, B and C are constants. Using these values in Eqs.(l7j)- (fTUl) . it 
follows that 

e 2r*(B-A)+2C( 2B _ ^ = g^ + gl ^ ^ 
e 2r\B-A)+2C^_ AB 2 r 2 + ^^2 + ^ = g^ _ ^ ^ 
e 2rHB-A) + 2C {4B 2 r 2 + ^ = g ffft + ^ ^ 
c 2r 2 (B-A)+2C 

a = - {re r2{A - B ^ c E)' . (18) 

Now we impose the central and boundary conditions on E(r) and p r (r) re- 
spectively as follows: 

E(0)=0, Pr (a) = 0, (19) 

where a is a positive constant and r = a is the interface of the charged fluid 
and vacuum (i.e., boundary of the cylinder). We apply central conditions to 
the system of Eqs. (fl5]) - (fT7|) . it follows that 



This shows that p r (r) = Pt{r) at r = 0, hence the anisotropy of the cylinder 
vanishes at the center. Applying the boundary conditions to Eqs.( lT5|) -( lT7)) . 

we get 

e 2a ^ B - A ^ +2C (-2B 2 a 2 + B + 2ABa 2 ) , s 

P{a) = ^ , (21) 

e 2a ^ B - A ^ 2C '(A + ABa 2 ) 
Pt(a) = ^ L , (22) 

E 2 {a) = e 2a2{B - A)+2C {AABa 2 + 2A - AB 2 a 2 ). (23) 
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The general expressions for p t and E 2 from Eqs. fll5l) - fll7l) are 

Pt(r) = - + P, (24) 

E 2 {r) = e 2r2{B - A)+2C (2B-2A)-87rp. (25) 



3 Models for Equations of State 

The general form of EoS is 

p r = p r (p,a 1 ,a 2 ), (26) 
where a\ and a 2 are parameters constrained by 

Pr (0) = p r [p(0) , a x , a 2 ) , = p r [p{a) , a x , a 2 ) . (27) 
Adding Eqs.f lTo"j) and f[T6"j) . we obtain 

e 2r 2 (B-A)+2C 

p + p r = (-4B V + 2B + AABr 2 ) = Z(r). (28) 

87T 

This equation may be used with the assumed EoS to find p and p r . The 
corresponding value of p will be used in Eqs. ( 1241) and ( 1251) to evaluate ^ and 
E 2 respectively. In the following we discuss three types of EoS: 

3.1 The Linear EoS 

This is given by 

p r = «i + a 2 p, (29) 

where cti, a 2 are constants. Using this EoS, we obtain expressions for 
p, p r , p t and E 2 as follows: 

e 2r*(B-A)+2C,_ AB 2 r 2 + 2B + AABr 2}_ 87rai 

87r(l + a 2 ) 

Q , 2e 2r 2 (B-A)+2C(_ 4B 2 r 2 + 2j B + 4ABr 2 ) + 871^ , . 

Pr = — • 31 

1 + OL 2 
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Using Eq.( 13"Uj) in Eqs. (l2"l|) and fl25]) successively, we get 

e 2rHB-A)+2C^ 2B _ AA _ 45V 2 ) (1 + a 2 ) + AB 2 r 2 - AABr 2 - 2B) 

8tt(1 + a 2 ) 

(32) 



Pt 



1 + a 2 



E 2 



e 



2r 2 (B-A)+2C 



{(2B - 2A)(1 + a 2 ) + 4S 2 r 2 - AABr 2 - 2B) 87r«i 



1 + a 2 1 + a 2 

(33) 



The values of constants «i and a 2 are found by solving Eqs. (|27|) and ( 129|) as 

P( a )Pr(0) Pr(0) 
1 " P(0) - P(a) ' " 2 p(0) - p(a) • ^ } 

3.2 The Nonlinear EoS 

The nonlinear EoS is given by 

Pr = ft + % (35) 

where n 7^ — 1 and are constants. It is a modification of the Chaplygin 

gas EoS used by Bertolami and Paramos [2 1 J to describe neutral dark stars. 
For n = 1, Eqs.fl28D and @5) lead to 



p = /(r)-A± ^(/(r)-^) 2 -^ (3g) 
Substituting Eq. (15Bj) in Eqs. (l3"5"j) . ( 12"4|) and ([2"5|) respectively it follows that 

Pr = Pl+ 2(32 =, (37) 

Kr)-p 1 ± y /(l(r)-p 1 )*-4p 2 ' 

e 2r*(B-A)+2C( B + 2£ >2 r 2 + 34) 
ft = 47T 

£ 2 (r) = e 2r2(B ~ A)+2C (2S-2A) -4yr(/(r) -A 

± vfOFl), (39) 



where l(r) is given by Eq. fl28j) . The constants fii and fa are found from 
Eqs. ([33]) and (J2ZD as 



p(0)-p(a) p(0)-p(a) 1 

Equation (I36I) implies that /3 2 must be negative so that each root in this 
equation has definite sign which will correspond to a positive definite matter 
density. 



3.3 The Modified Chayplygin Gas EoS 

This EOS has the following form 

Pr = JlP+J, (41) 

where 71, 72 are constants. This is used to describe static, neutral, phantom- 
like sources |22j. Using this EoS with Eq.f l28j) . we get 



= /(Q + y^- 4(1 + 71)72 
9 2(1 +71) ' 1 ] 

Substituting this value of p in Eq.( l4"T|) as well as in Eqs.( l24"j) and ( 125]) suc- 
cessively, it follows that 



Pr = 7i 



'l {r)±^l{rf- 4(1 +71)7 2' 
2(1 + 7i) 



Pt 



+ ™ 1 + *> (43) 

l(r)± v//(r) 2 -4(l + 7l ) 72 

e 2r\B-A)+2C(_ B + 2B 2 r 2 + 2A) 

Air 



l(r)± y//(r)2- 4(1 + 71)72 , . 

2(1 + 71) ' 1 ] 

E 2 = e 2r 2 (B-A)+2C( 2B _ 2A) 

4 . ( MWggES^y (45) 
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where 71 and 72 are 



/9(0)p,(0) p(0) Pr (0)p(a) 2 
71 p(0P-p(a)2' 72 - p(oy-p( a y { ) 

For (1 + 71) < and 72 > or (1 + 71) > and 72 < 0, Eq.flfS]) yields roots 
of definite sign which will again correspond to the positive matter densities. 
It is clear that the positive or negative matter densities depend upon the 
positivity or negativity of the constants fii, f3 2 , 71 and 72, i.e., the roots of 
Eqs.( |36|) and (1421 respectively. 



4 Matching Conditions and Adimensional Mat- 
ter Sources 

Here we take the charged static cylindrically symmetric spacetime as an 
exterior region given by 



ds 2 = N(r)dt 2 - j^dr 2 - r 2 d6 2 - r 2 d^ 2 , N(r) = ^ - (47) 

where q and m are charge and mass respectively. Using the transformation, 
d^) = J- — d(j), this takes the form 

ds 2 = N(r)dt 2 - -r^--dr 2 - r 2 d6 2 - -r^-dcj) 2 . (48) 
v ; N(r) N(r) v ; 

2 2 

With the radial transformation r = m ~ q , it becomes 

q 2 2m i 1 2mq 2 



im z — q A \2 1 2m ,i _i_ 2mq z v T ' v ' 

^ 172 ) Z12 "T73" r /3 



To match the interior metric ([2]) with the exterior (149]) . we impose the 
continuity of g 00 , 9i\ an d ^f^ 1 across a surface at r' = a by using the procedure 
j. In our case, this yields the following expressions for A, B and C in terms 
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of adimensional parameters rj = — and x = as 

_ ln^ 2 - x 2 ) , v 

" a 2 (?7 2 -x 2 ) 2 ' ( j 

1 1 V 2 
R = f I * ^ 

a 2 {r] 2 - X 2 ) 2 A{ X 2 -2rf + r 1 x 2 ) 

C = - + — -= ^ + In -^=JL ~ X (52) 

2 4( X 2 - 2rf + VX 2 ) yfi? - Irf + VX 2 

We can see from Eq. (fT4l) that A and .B have dimension of length~ 2 and C 
is dimensionless. It is very important that the field equations can eventually 
be expressed in terms of these adimensional constants and the dimension- 
less radial coordinate x = -. Here adimesionality is denoted by hats. We 
assume that the interior of the fluid cylinder is described by x G [0, 1). 
We reformulate all models as adimensional models where we are denoting 
A = a 2 A, B = a 2 B, p = a 2 p, p r = a 2 p r , p t = a 2 p t , E 2 = a 2 E 2 and a = a 2 a. 
The quantities which are originally dimensionless are denoted by the actual 
symbol. The central and the boundary conditions at x £ [0, 1) become 

(B-A)e 2C , N 

m = { J > (53) 

Ae 2C 

ft-(O) = (54) 

Ap 2C 

MO) = (55) 

e 2 ^-^ +2C (2S 2 — B — 2AB) , . 

P(l) = 4^ K (56) 

Ml) = J + ] , (57) 

E 2 {1) = e 2 ^~ A)+2C {AAB + 2A-AB 2 ). (58) 

The Adimensional linear EoS model 

The adimensional linear EoS is given by 

p r = &i + a 2 f>. (59) 
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The corresponding quantities will become 

e 2x2 ^- A ^ +2C (-AB 2 x 2 + 2B + AABx 2 ) - 87^ 
P ~ 8tt(1 + a 2 ) ' 

87^ + a 2 e 2x2 ^' A ^ +2C (-AB 2 x 2 + 2B + AABx 2 ) 



Pr = 



l + a2 



where 



2x 2 (B-A)+2C 

Pt = - (-2B + AB 2 x 2 + AA) + p, 

E 2 = e 2x2 ^-^ +2C (2B-2A)-8np, 



p(i)p(Q) Mo) 

OL\ = —-T7^ T7TT, «2 ~~ 



p(0)-/3(l)' Z p(0)-p(l)" 

The Adimensional Nonlinear EoS Model 

Here we have 

k M^ 2 
Pr = Pl + — • 

P 

For this model, the corresponding quantities take the form 

i(x) -$i± J{i{x)-h) 2 -±h 
p = 2 ' 

* = 4^ L+p > 

E 2 = e 2x2 ^- A ^ 2C (2B-2A)-8irp, 
where (3i and (3 2 are 

h p(0)p r (0) * P(0)p(l)p r (0) 

Pl p(o)- p (i)' P2 " p(o)- p (i) • 

The Adimensional Modified Chaplygin EoS Model 

The adimensional modified Chaplygin EoS model 

- . 72 

Pr = 7lP + — 
P 
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yield the following quantities 

l{x) ± y / f(x) 2 -4(l + 7l )7 2 



2(1 +71) 



(71) 



where 



e 2x 2 (B-A)+2Cfg _^ 2B 2 x 2 + 2A) , , 

ft = 4^ J -+P, (72) 

E 2 = e 2x2 ^~ A > +2C (2B-2A)-8irp, (73) 

P(0)Pr(0) « p(0)p r (0)p 2 (l) 

71 p2(0)-/)2(l)' 72 " ^(O)-^(l)' 1 J 

The adimensional proper charge density for all the above three models will 
become 

p 2x 2 (B-A)+2C J „ . „ 

a = ±( xe * 2 (B-A) + 2Czy (?5) 

2ttx dx y ' y ' 

5 Some Features of the Models 

In this section, we discuss some insights of the three models. The exterior 

2 

metric p7|) implies that singularity occurs at r = 0, It is regular ev- 

2 

erywhere except at r = [23J. The surface with r = ^— describes a right 
singular circular cylinder. Using Eqs. (!50|) - (l52|) in p(0) of Eq. (l20|) . it turns out 
to be positive. Following [T2], we have t] = m 1.147 for M to be the solar 
mass and a = 1.48/cm. For this value of r], p(0), p(l), pt(0), Pt(l), -E'(l) 
and p r (0) turn out to be functions of x only. Similarly, the expressions for 
&i, a 2 , (3 2 , 7i and 72 (EoS parameters in adimensional version) also 
depend only on \- The analysis of A, B and C implies that the values of x 
are restricted by the values of 77 such that x < 77 for 77 > 0. 

Figures 1 and 2 indicate that the central density p(0) and central pressure 
p r (0) are monotonically decreasing with increasing values of x f° r V — 1-147 
and x £ (0, 0.7]. We see that p r (0) > only when x £ [0, 0.57), hence we take 
the maximum charge x = 0-56 to discuss our models. Figures 3-5 show that 
$2 < 0, 1 + 71 > and 72 < for the same value of 77. We can have positive 
definite roots from Eqs. (I66j) and (j7TJ . hence we can analyse the models with 
positive matter density p(x). It is obvious that x = implies no charge. 
If we increase x, it increases repulsive electrostatic forces and consequently 
pressure and density change. 
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Now we explore the behavior of matter sources for our models at different 
charges. Firstly, we investigate the sources with linear EoS. Figures 6-9 
show the corresponding matter density, radial pressure, tangential pressure 
and electric field intensity. The graphs of density and pressure show the 
increasing behavior while the electric field intensity is decreasing with the 
increasing values of x- Further, E(x) decreases at every point in the interval 
x G (0, 0.1] with increasing %■ Figure 10 shows the behavior of charge density 
a which is unbounded for each value of x- 

The analysis for models with nonlinear and Chaplygin gas EoS indicates 
that these models are similar to the model satisfying linear EoS. These models 
correspond to the decreasing matter densities and pressures. Each EoS affects 
the dependance of the measure of anisotropy 5 on x. The only difference 
arising from the three models is about the measure of anisotropy 5 = pt — p r . 
Figure 11 displays the anisotropic parameter for the model corresponding 
to the linear EoS which is increasing with the increasing x- Figures 12 
and 13 show the anisotropic parameters corresponding to the nonlinear and 
Chaplygin gas EoS respectively. From figures, we see that the anisotropic 
parameter for nonlinear EoS is increasing while for the Chaplygin gas EoS, it 
is decreasing with increasing x- I n addition, Eq. (fl8|) contains proper charge 
density a which implies that 

rae 2 ^ 2v dr = (76) 
where 

q {r) = Anre u -^ / rae 2 ^ 2u dr (77) 
Jo 

which is the net charge inside the cylinder of radius r. 

For our cylindrically symmetric fluid models, Eqs.fl7])-f lTU]) represent the 
basic source parameters. We formulate table I by computing adimensional 
values for these basic sources. Approximate numerical values for central 
density p(0), central pressure p(0), tangential pressure PtiX) and electric field 
intensity E 2 (l), are shown in table I. We note that the maximum value of 
p(0) and the value of p r (0) are smaller than maximum of j3 t (l). These max- 
ima correspond to zero net charge sources. The source variables p(0), ]3 r (0) 
with maximum charge (x = 0.56) have been changed while pt{0) changes its 
sign at the maximum charge. To compensate the stronger electric repulsion, 
this sign inversion of tangential pressure is necessary. This is summarized in 
the following table: 
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Table I. Approximate numerical values of some quantities. 



X 


P(0) 


p(0) 




E\l) 





0.0221364 


0.00651431 


0.00830257 


0.266397 


0.56 


0.0172785 


0.0000738598 


-0.0000726 


0.370451 



Now we discuss some consequences of our results. 
5.1 Stability 

Bertolami and Paramos [21] argued that if the generalized Chaplygin gas 
tends to a smooth distribution over space then most density perturbations 
tend to be flattened within a time scale related to their initial size and the 
characteristic speed of sound. One of the important physical acceptability 
conditions for anisotropic matter is that the squares of radial and tangential 
sound speeds (u 2 = and u\ = -j^) should be less than the speed of light 
|25j . We explore it for the linear EoS model. The graphs of the squares 
of radial and tangential sound velocities are shown in Figures 14 and 15 
respectively. These indicate that v% is independent of x and decreases with 
increasing x while uf monotonically increases with increasing x and also 
increases with increasing x for fixed x. For three particular values considered 
here, these parameters satisfy < v% < 1 and < v% < 1 everywhere within 
the charged fluid. 

Now we use Herrera [25] and Andreasson's [26] approach to identify po- 
tentially unstable or stable anisotropic matter configuration. According to 
their approach, \uf — uf \ < 1 as shown in Figure 16. This implies that 

(?) - 1 < u\ - u 2 T < 0, (ii) < u\ - u 2 r < 1. 

The first expression corresponds to the potentially stable model which is 
obvious from Figure 17, while the second expression corresponds to the po- 
tentially unstable model. We note that our model satisfy the potentially 
stable condition. If the graph of u\ — u 2 keeps the same sign everywhere 
within a matter distribution, there will be no cracking and the system is 
stable. If there is a change of sign then it is alternating potentially unstable 
to stable region within the matter distribution and vice versa. 
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5.2 Energy Conditions 

The energy conditions of the charged anisotropic fluid, the weak energy con- 
dition, the strong energy condition and the dominant energy condition are 
satisfied if and only if the following inequalities hold: 



P+Pr 


>o, 


(78) 


E 2 


>o, 


(79) 


E 2 

P + Pt + ^ 


>o, 


(80) 


E 2 

p + p r + 2p t + — 


>o, 


(81) 


E 2 E 2 

P+2 \Pr ~ 


>o, 


(82) 


E 2 E 2 

P + Q Pt + ^ 

8tt Stt 


> o. 


(83) 



Figures 18-23 indicate that these inequalities hold for each x G [0,1]. 
5.3 The van der Waals (VDW) EoS 

Lobo [19] introduced this type of bounded source for the sake of cosmology. 
Here we express this approach with VDW EoS given by 

Pr = -^-u 2 p 2 . (84) 

1 - LU 3 p 

This is used to describe dark matter and dark energy as a single fluid. It 
is assumed that the interior and exterior metrics are joined and E(0) = 
= p r (a). Further, p r is found at the center and boundary of the charged 
cylinder. Using Eqs. (l84"]) and (I3"0l) . it follows that 

uj 2 u 3 p 3 - (co 2 + u 3 )p 2 + [1 + ui + uj 3 l(r)]p - l(r) = 0, (85) 

where U\ and u 2 are functions of p(0), p(a), p r (0) and u 3 . We do not expect 
any interesting consequences from this equation, hence we leave it here. 
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5.4 Equilibrium Condition 

Now we discuss the variation of the net charge corresponding to different 
forces compatible with equilibrium configuration for our models. In partic- 
ular, when pressure gradients tend to zero and the charged fluid is more 
diluted, then what is the behavior of gravitational and other forces. Using 
Tolman-Oppenheimer-Volkov equation, we obtain 

(p+p r )(e^(p r -i? 2 )-e 4 - 2 ^)-^ + ^^ + r(e 2 ^)fe-p r ) = 0. (86) 

r ar r z 

This provides the equilibrium condition for the charged fluid elements subject 
to different forces. Here q = q(r) as given in Eq. flTTj) . In adimensional version, 
we can write 

F 1 + F 2 + F 3 + F 4 = 0, (87) 



where 



2Bx 2 +2C/~ fr2\ S ABx 2 -2Ax 2 +AC\ 



F l = (p + p r )(e 2Bx+2C (p r -E 2 )-2Axe 

* ■ ~% < 89 » 

F 3 = aEe x2(B ~ A+c) , (90) 
F 4 = xe 2x2 ^- A+2C \p t -p r ). (91) 

The graphs of these forces with linear EoS at x — and x — 0.56 are shown in 
Figures 24 and 25 respectively which indicate that the charge has a negligible 
effect on these forces, hence we obtain a static equilibrium. For x — 0, F 4 
point outwards at every x G (0, 1] and F 2 is along the x-axis. The electric 
force F3 acting on the fluid elements with unbounded a located at x — 1 is 
infinite. This is the weakest force because it changes sign at x w 0.38 and also 
the force F\ changes sign at x ~ 0.9 which is due to gravity. When x — 0.56, 
the electric force is still unbounded and infinite. This unboundedness of the 
force F 3 and sign inversion of the force iq is essential for the configuration 
of our static, charged anisotropic model with linear EoS. 

6 Outlook 

The main purpose of this paper is to investigate the solutions of the cou- 
pled Einstein-Maxwell field equations for the static cylindrically symmetric 
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spacetime. For this purpose, we have used charged anisotropic fluid with EoS 
in the light of Victor et al. [12] procedure developed for static spherically 
symmetric spacetime. In particular, the linear, nonlinear and Chaplygin EoS 
have been used. It is mentioned here that the linear EoS corresponds to the 
electrically charged isotropic strange quark stars on the basis of MIT bag 
model [27]. Our model with linear EoS also corresponds to the model by 
Victor et al. [12] for the charged anisotropic spherically symmetric fluid with 
the same EoS. We know that the nonlinear and Chaplygin EoS are used to 
describe non-static neutral gravitational isotropic fluids. Here we are taking 
static charged anisotropic fluid as we would like to explore the interior re- 
gions and the fluid vacuum interfaces of the charged anisotropic cylindrically 
symmetric stars with these EoS. 

We have used the assumptions of Karori and Barua to explore the charged 
anisotropic static cylinder. Our models with nonlinear and Chaplygin EoS 
correspond to the dark matter and dark energy with constant matter den- 
sities and pressures [28]. Delgaty and Lake [29] proposed some physical 
conditions acceptable for perfect fluids, i.e., regularity of the charge at the 
origin, positive matter density and pressure, decreasing matter density and 
pressure with increasing r, causal sound propagation and smooth matching 
of internal and external metrics at boundary of the source. Burke and Hobill 
[30] added one more condition that sound velocity must be monotonically 
decreasing with increasing r which was imposed on spherical perfect fluid 
[3TJ [32]. We have found that our cylindrical models satisfy most of the above 
physical conditions. In our case, conflict may arise for increasing tangential 
sound velocity, constant radial velocity and negative tangential pressure for 
higher values of x- 

Finally, we would like point out that we have also explained charged 
anisotropic static cylindrically symmetric models in our recent paper [T7] by 
using Thirukkanesh and Maharaj approach. In that paper, we have found 
that the charge distribution as well as E become singular at r = 0. However, 
here we have found non-singular behavior of these quantities. 
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Figure 1: Variation of central density p(0) with x- 
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Figure 2: Variation of central pressure p r (0) with x- 



-0.00002 
-0 .00004 
-0 .00006 
-0 .00008 
-0 . 0001 
-0.00012 



0.1 0.2 0.3 0.4 0.5 



Figure 3: Variation of parameter /3 2 with %■ 
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Figure 4: Variation of parameter 71 with x- 




Figure 5: Variation of parameter 72 with \. 
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Figure 6: Variation of matter density p(x) with x for different values of x- 
Red, green and purple curves correspond to x — 0, 0.3, 0.56 respectively. 
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Figure 7: Variation of radial pressure p r (x) with x for different values of \. 
Red, green and purple curves correspond to x — 0, 0.3, 0.56 respectively. 




Figure 8: Variation of tangential pressure pt(x) with x for different values of 
X- Red, green and purple curves correspond to x = 0, 0.3, 0.56 respectively. 




Figure 9: Variation of Electric field intensity E(x) with x for different values 
of x- Red, green and purple lines correspond to x = 0, 0.3, 0.56 respectively. 
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Figure 10: The charge density a(x) for different values of %. Red, green and 
purple curves correspond to x — 0, 0.3, 0.56 respectively. 
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Figure 11: Measure of anisotropy 5 for the first model with x for different 
values of x- Red, green and purple curves correspond to x — 0, 0.3, 0.56 
respectively. 
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Figure 12: Measure of anisotropy 8 for the second model with x for differ- 
ent values of x- Red, green and purple lines correspond to x — 0, 0.3, 0.56 
respectively. 
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Figure 13: Measure of anisotropy 5 for the third model with x for differ- 
ent values of \. Red, green and purple lines correspond to x — 0,0.3,0.56 
respectively. 
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Figure 14: Variation of radial sound velocity u r (x) with x for different values 
of x- Red, green and purple lines correspond to x = 0, 0.3, 0.56 respectively. 
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Figure 15: Variation of tangential sound velocity ut(x) with x for different 
values of \- Red, green and purple curves correspond to x — 0, 0.3, 0.56 
respectively. 
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Figure 18: Weak Energy Condition 1 
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Figure 19: Weak Energy Condition 2 
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Figure 20: Strong Energy Condition 1 
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Figure 21: Strong Energy Condition 2 




Figure 22: Dominant Energy Condition 1 




Figure 23: Dominant Energy Condition 2 
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Figure 24: Four different forces Fi(x), F 2 (x), F 3 (x) and F±{x) corresponding 
to red, green, purple and blue curves respectively are acting on the fluid 
elements in static equilibrium for x = 0. 




Figure 25: Four different forces Fi(x), F 2 (x), F 3 (x) and F A (x) corresponding 
to red, green, purple and blue curves respectively are acting on the fluid 
elements in static equilibrium for x — 0.56. 
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